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Abstract 

The eigenvalue problem of the Hamiltonian of an electron confined to a plane 
and subjected to a perpendicular time- independent magnetic field which is 
the sum of a homogeneous field and an additional field contributed by a 
singular flux tube, i. e. of zero width, is investigated. Since both a direct 
approach based on distribution-valued operators and a limit process starting 
from a non-singular flux tube, i. e. of finite size, fail, an alternative method 
is applied leading to consistent results. An essential feature is quantum me- 
chanical supersymmetry at g = 2 which imposes, by proper representation, 
the correct choice of "boundary conditions". The corresponding representa- 
tion of the Hilbert space in coordinate space differs from the usual space of 
square-integrable 2-spinors, entailing other unusual properties. The analysis 
is extended to g 7^ 2 so that supersymmetry is explicitly broken. Finally, 
the singular Aharonov-Bohm system with the same amount of singular flux is 
analysed by making use of the fact that the Hilbert space must be the same. 



1 Introduction 



The parameter space of the system that is considered interpolates between two very 
different systems of fundamental importance. 

On the one hand, the motion of an electron in a plane perpendicular to a homo- 
geneous magnetic field [U@ has found an important application in quantum Hall 
physics [rj|[Q||§]. Its most striking feature is the vast degeneracy within the discrete 
Landau levels that makes quantum Hall phenomena possible. It is of interest how 
this system is altered by inhomogeneities of the magnetic field. An additional sin- 
gular flux tube appears to be a minimal modification of the homogeneous field ||, 
while in general inhomogeneous magnetic fields are notoriously difficult to handle. 

On the other hand, since the work of Aharonov and Bohm |7| the physics of a 
magnetic flux tube alone has been investigated under a variety of aspects. Con- 
sidering particles of spin 1/2 M, the system is of interest as a simple application 
of mathematical index theorems |T(| |L1] , interrelated with anomalies in quantum 
field theories [I2|[I3|[I3| and the fact that the Pauli equation exhibits a quantum 
mechanical supersymmetry for g = 2 [[OJ. Furthermore, it is particles of spin 1/2 



that problems occur for, when we go over from the non-singular, i. e. flux tube of 
finite width, to the singular case, i. e. flux tube of zero width. While for the non- 
singular case the standard quantum mechanics works well, the singular case raises 
the question of the correct boundary conditions at the location of the flux tube. 
Standard boundary conditions, e. g. of the usual Dirichlet type, are not sufficient 
to characterize the behavior of the eigenfunctions of the Hamiltonian ||16| ||17|| . This 
is due to the fact that functions of different boundary behavior are mapped into 
each other by the supersymmetry, which is thus not represented automatically by 
respecting simple boundary conditions. This difficulty makes one suspect that the 
transition from the non-singular to the singular flux tube might not be continuous. 
This point among other questions is obscured by the fact that a flux tube alone, 
yielding the free particle behavior as we go to infinity, leads to a continuum in the 
spectrum of the Hamiltonian entailing non-normalizable eigenfunctions. 

We supply the singular flux tube of flux a$, where a is real and $ = hc/\e\, 
with an additional homogeneous magnetic field and solve the eigenvalue problem of 
the Hamiltonian. Technically our procedure deviates from standard methods con- 
siderably. Some features that we are familiar with from other quantum mechanical 
systems are not found, while other features that are unusual in a quantum mechan- 
ical context do occur. Therefore, we want to clarify some issues relevant to our 
system and point out some peculiarities of the system in advance. 

The integral part of a is not a gauge effect. The similarity to a gauge trans- 
formation merely shows that an inhomogeneity constrained to a single point leaves 
many features of the system intact. Yet, systems differing by integral quanta of 
singular flux are different. 

We solve the non-relativistic Pauli equation and not the corresponding Dirac 
equation. A relativistic treatment would introduce the Compton wave length Ac 
as an additional scale to the system. An infinitely thin flux tube would run into 
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trouble with the consistency of one-particle quantum mechanics by introducing a 
scale falling below Ac. While a magnetic field does not accelerate a charged particle, 
a gradient exerts a force on the magnetic dipole. We could be confident about the 
one-particle description if the radius of the flux tube were at least of the magnitude 
M 1 ^c- Otherwise, pair production due to the inhomogeneity of the magnetic 
field had to be taken into account. Staying non-relativistic we avoid any conflicts. 
Clearly, this means that a singular flux tube idealizes a real flux tube of minimal 
radius |a| 1//2 Ac- 

The quantum system will not be realized in coordinate space before the limit 
R — > is taken, where R denotes the radius of the tube to which the flux a$ 
is confined. In the regularized coordinate representation we quantize the system 
and calculate matrix elements. The regularization is an intrinsic necessity in order 
to handle unavoidable products of singular objects at the flux tube consistently, 
since matrix elements contain products of both singular operators and singular wave 
functions. Supersymmetric pairing turns out to be possible for eigenfunctions only 
due to an arbitrariness following from the regularization. Also, this regularization 
makes possible a deviation of the Hilbert space from the space of square-integrable 
2-spinors L2 <8> C 2 for \a\ > 1. In particular, singular supersinglet states occur with 
properties similar to those of a classical point particle in two dimensions. 

Usually, after solving an eigenvalue problem of a Hamiltonian we verify that the 
set of eigenfunctions supplies a resolution of unity in the space of square-integrable 



functions [18]. For our singular system, however, we do not know a priori what 
the Hilbert space should be, i. e. what space of functions we should find a resolu- 
tion of unity for. Conversely, we use the orthonormal set of eigenfunctions of the 
Hamiltonian to span the Hilbert space. The completeness of the eigenfunctions of 
the Hamiltonian within the Hilbert space is an indispensable requirement that is 
accounted for in this way by the very construction. 

As opposed to non-singular systems, in the presence of a singular flux tube the 
representation of supersymmetry is not automatic by a standard treatment using 
local boundary conditions |13|] . Instead, insisting on a preservation of supersymme- 
try uniquely determines the correct eigenfunctions of the Hamiltonian and, hence, 
the correct Hilbert space. 

We specifically treat the electron with charge — |e|. The considerations apply to 
any elementary, charged, massive particle of spin 1/2 by obvious modifications. 



2 Supersymmetry 



The following statements JTI| |20] [[Uj rely on the assumption that all quantities are 
sufficiently differentiable and non-singular so that their products are defined. 

We consider the Pauli Hamiltonian in the x-y-pleaie with p z = 0, g = 2, and a 
perpendicular magnetic field B = B z (r)e z , where r — (x,y) 

H = ^-(v 2 x + v 2 y )+ l -^B z S z , (1) 
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with the velocities v x = (p x + ^-A x )/M and v y = (p y + ^A y )/M. Angular momenta 
are always given in units of Ti. We define the superoperators Q = S + V and = 
with V = (M/2) 1 / 2 ^ + iv x ), = {M/2fl 2 {v y - iv x ) and S± := S x ± iS y . 
The Hamiltonian and the superoperators obey the supersymmetry algebra 

Q 2 = (g t ) 2 = 0, H = QQ^ + Q^Q. 

^[H,Q] = [H,Q^] = 0. (2) 

In addition we have 

[S X ,Q]=Q, [S*,Qi] = -Q\ [S z ,H]=0. (3) 

Hermitian superoperators are Q = Q + = (2M) l l 2 (S x v y — S y v x ) and = 
i{Q ] -Q) = {2M) l l 2 {S x v x + S y v y ). 

We can immediately draw important conclusions from supersymmetry. On the 
one hand, the measured energy is 

mH\^) = (Qtf|Qtf) + (Qttf|Qttf) > . (4) 

On the other hand, any eigenstate of the Hamiltonian with eigenvalue E > is 
doubly degenerate. By (|3|) we can always arrange the partner states to be eigenstates 
of ^ such that 

Q\E,a = -1/2) ~ \E,a = +1/2), 

Q^\E, a = +1/2) ~ \E, a = -1/2), (5) 

where S z \a) = a\a) with o = ±1/2, if E > 0. The eigenstates of H with E = 0, 
being annihilated by the superoperators 

Q\E = 0) = Qi\E = 0) = 0, (6) 

are supersinglets. 

3 Insufficient approaches to the problem 

We give a brief sketch of two obvious approaches in order to work out their short- 
comings resolved by the correct method in the next section. 

3.1 The direct approach 

We consider the magnetic field in the form B(r) = (B + a$<5 2 (f))e 2 ||. B > 
is constant everywhere. The corresponding vector potential is A(r) = (Br/ 2 + 
a$/27rr)e ¥ ,, where r and (p denote polar coordinates in the x-y-pla.ne] i. e. the 
flux tube part is a pure gauge for r ^ locally. We measure lengths in terms of 
A = {$>/nB) l l 2 putting f = r/A. Energy is measured in units of Tluj with the Larmor 
frequency to = \e\B/Mc of the homogeneous field. The superoperators are measured 
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in units of (hcu) 1 ^ 2 . For a = we obtain the operators for the homogeneous magnetic 
field. 

The Hamiltonian is 

1 9 2 \ 1 „ 2 2 9 a za; 9 a 2 r/~\n / 7 \ 

Q = ~4 ^97 r 97 + ^9^J + 4 r ~29^ + 2+ 2 _ 2^9^ + 4^ + 2f (r) " ( } 

Since [# a ,Lj = [if a ,Sy = 0, we put ^ E ,a,m{f, <p) = ipE,a,m(f)e im,fi Ca in the eigen- 
value equation, where m is integer and ( a are unit eigenspinors of S z . For f 7^ the 
eigenvalue equation is the same as for a homogeneous magnetic field by the substi- 
tution m — > m + a. From the eigenf unctions of the homogeneous field a = 0, we 
know the solutions for a 7^ in the entire plane except at f = 0. The delta function 
allows only eigenfunctions to be continued to f = that vanish there. Accordingly, 
the solutions are 



*E,<7,m(r, (p) = 



(E - a - 1/2 - (m + a)9(m + a))\ 
T(E - <x + 1/2 - (m + a)6(-m - a)) 

a-l/2-(m+a)e(m+a) ' 



v fl m +«l j\ m + a \ (r 2 \ P ^r ■p imi P/- 



with E — o — 1/2 — (rn + a)6(m + a) = 0, 1, 2, ... and |m + a\ 7^ 0. For consistency 
with supersymmetry it is necessary that each eigenfunction with E > find a 
superpartner. The superoperators 

and 

«t=S-ltf = S-£( »-'» +r - + HV (10) 



however, raise or lower the leading power for small f by one unit such that an 
eigenfunction (§) that is 0{f u ) with < u < 1 might be mapped to a superpartner 
that is 0(f u_1 ), which is not an admissible eigenfunction of H a , being in conflict with 
the delta function. Insisting on supersymmetry, we have to exclude all eigenfunctions 
in (H) of which the superpartner is not contained among the eigenfunctions in (H) 
as well. 

For each spin component for non-integer a one value of (m + a) is missing. For 
integer a^O even two values of (m + a) are missing. The set of eigenfunctions of 
the Hamiltonian certainly does not provide a complete basis for the space of square- 
integrable 2-spinors L2 ® (D 2 if a 7^ 0, because of the missing m- values. Although it 
appears strange at first sight, this is no reason to reject this approach. One should 
also recall that the canonical orbital angular momentum L z is not a gauge-covariant 
operator and, therefore, is not observable. The eigenstates of this direct approach 
are contained among the correct eigenstates displayed in figure 1. However, for 
a > the column of states left to the defect line and for a < the column right to 
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the defect line is missing. Also the singular states given in the figure are not visible 
by the direct approach. Incidentally, if there were no pairing, this direct approach 
would still yield a vacancy for a = ±1, ±2, ... . 

It is not obvious that {Q a , Q^} = H a , i. e. that the Hamiltonian with the delta 
function as in is indeed a result of the anti-commutator. This result can be 
obtained using complex coordinates Z = fe tip , defined also for f = by using the 
identity 

~=^m, (id 



valid on functions that are continuous at f = |2T| . The identity is derived by using 
a regularization that cannot be avoided. In general, products of singular operators 
such as Q a , Q^, and H a require a regularization. 

Distribution-valued operators require sufficiently well-behaved wave functions 
playing the role of test-functions so that matrix elements are always defined. This 
is, implicitly, the idea behind the above approach to the eigenvalue problem of 
H a , where the possibility that the wave functions themselves could exhibit singu- 
lar behavior is excluded by inconsistency with the formulation of the problem in 
terms of distribution-valued operators. An extension beyond this setting requires 
a regularization for both the operators and the wave functions in order to make 
the corresponding products defined. A treatment of the eigenvalue problem of H a 
based on such a regularization will have to reproduce all eigenfunctions of the above 
method, which are certainly correct. It could still, however, reveal other more sin- 
gular eigenfunctions. 

3.2 The singular flux tube as a limiting case of a flux tube 
of finite size 

Now we consider the system of a cylindrical flux tube of finite radius R. Apart from 
the homogeneous field throughout the plane, we add another homogeneous field for 
r < R contributing the flux oe<&. The magnetic field is 

B(r) = [B + 0(R - r)a$/iTR 2 }e z , (12) 

where 6{R — r) = 1 for r < R and otherwise. The corresponding vector potential 
is 



1(f) 



Br „, , a$ „. . «$r 
R) + 9(R-r] 



2 v J 2rrr y ' 2kR? 



(13) 



We employ r = r/X and R = R/X with the magnetic length A of the homogeneous 
field outside of the flux tube and the corresponding tkv as the unit of energy. The 
Hamiltonian is given by 

H -- (-—f— ——\ + -f 2 - -— + S 

4 \f df df f 2 dip 2 J 4 2 dp 2 



5 



+ 6>(r - i?) - — — I + 9{R - r)—[ -^r 2 + -r 2 - -— + S z 

and the superoperators are given by 



(14) 



Qc 



s+v n = s 4 



-lip 



d id a ~ / r 1 

ttz - + r + - + ol6{R - r) U- - - 
or r r \R 2 r 



Ql 



a 
f 



/ r 




Vl2 


r / 



(15) 



9 z 9 

The supersymmetry algebra @ is fulfilled. By [i7 a , S z ) = [H a , L z ) = we put 
^E,o-,m(^, V 9 ) = ipE^mi^e' 1 ™ 11 ?^ for the eigenfunction. The eigenvalue equation of 
H a supplies differential equations of Kummer type both inside and outside of the 
flux tube. The inside solution has to be regular at f = 0, while the outside solution 
must decay as f — > oo. We match the two solutions by demanding continuity at 
f — R yielding 



E,a,m 



+9(R-r 



6(f - R)f m+a e—U(m + a + a + 1/2 - E,m + a +l,r 



R m + a -\m\U( m + a + a + 1/2- E,m + a + l, R 2 )e h - 



x F x (mQ{m) + a + 1/2 - p , 2 
xe- (1+a/ ^ 2)f2/2 f |m| 1 F 1 (m^(m) + a + 1/2 



R 2 E 



R 2 + a 



1{2E ,l + \m\,R 2 + a 
1 + \m\, II 



a 
~R? 



r 2 ) 



tmip 



(16) 

with the Tricomi function U [p2 |. Still we have to demand continuity of the first 
derivative at f = R. This entails the following equation determining the eigenvalue 
E for given values of a, m, a, and R. 

52 U'(m + a + a + 1/2 - E, m + a + 1, R 2 ) 



R 



U{m + a + a + 1/2 - E, m + a + 1, R 2 ) 

fm#(-m) + a = 0. (17) 



xFx'imdim) + a + 1/2 - Jjg-, 1 + |m|, R 2 + a 



l F 1 (m6(m) + a + 1/2 - 



B?E 
R 2 +a 



1 + ImLi? 2 + a) 



The prime indicates differentiation with respect to the third argument. 

Since we are interested in the limit R — > we consider ( |i~7[ ) for small R and find 
the condition 

(±m ± a + a — \m\)iFi(mO(m) + a + 1/2, 1 + \m\, a) 

- 2a l F l '(m9(m) + a + 1/2,1 + |m|,a) = (18) 

up to terms that vanish as R — > if R m+a U(m + a + a+ 1/2 — E, m + a + 1, J? 2 ) = 
0(i? ± ( m+a )). First of all, we observe that ( jT8|) is independent of E. But considering 
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the various cases yields a further simple result. For a = — 1/2 and m < on the 
one hand, and a = +1/2 and m > on the other hand, fllTf ) is identically fulfilled. 
Otherwise fll8f) is never fulfilled. Consequently, if one superpartner is in agreement 
with QTSf) , the other necessarily violates (|I8|). Thus, supersymmetric pairs are not 
admitted and only supersinglets can occur. Therefore, only eigenf unctions with 
E = are possible. They have the particularly simple form 

V E=0 CT= _ 1/2 m (f) ~ [0(f-£)f-^ a >e^+0OR-f ) J R- a et e - (1+a/ ^ 2)r " 2/2 f : - m ]e 4m ^C-i/9, 

(19) 

for all R and a; with m = 0,-1, —2, ... They can also be found by directly solving 
Q a ^E=o(r, ip) = e=q(? i = 0. These eigenfunctions are the only ones that 
survive the limit R — > 0. Since they can always be normalized so that their norm is 
finite in the limit R — > 0, all of them have to be considered as correct. 

For a > these singlets have quantum numbers for which no eigenfunctions 
are found by the direct approach. By contrast, this approach has not provided the 
majority of eigenfunctions found to be correct by the direct approach. Thus, there is 
a mutual incompatibility of the two approaches of this section. The correct method 
will have to reproduce the safe results of both alternatives. 

We point out that the failure of the approach of a shrinking flux tube is not a 
result of its rough implementation. If the magnetic field profile is made smooth by 
a modification within the annulus between R — 7] and R, the limit R — > will lead 
to the same results. 

4 A consistent method 

4.1 The eigenvalue problem of the Hamiltonian 

The direct approach in 3.1 suggests that a regularization should be applied. How- 
ever, the most obvious way to regularize the problem, by regarding the system of 
the singular flux tube as the limiting system of a sequence of systems of shrink- 
ing flux tubes in 3.2, fails. Although each member of the sequence is a consistent 
quantum mechanical system, a limiting system does not exist in agreement with the 
obviously correct solutions of the direct approach. Therefore, the idea is to replace 
the sequence of physical systems at finite R by a sequence of unphysical systems 
at finite R. Only the limiting system is physical. In other words, for finite R the 
mathematical structures cannot be interpreted in a quantum mechanical sense. A 
consistent quantum mechanical interpretation is possible only in the limit R — > 0. 
This is comparable to the use of the dimensional regularization in quantum field the- 
ory, where the system may not be interpreted in physical terms before the regulator 
is removed. 

Before establishing the formal framework, we anticipate the essential points that 
suffice for a solution of the eigenvalue problem. We use the operators of the finite 
flux tube as given in 3.2 as their regularized form. Any matrix element of an operator 
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f2 is to be calculated as 

(E|fi|*) = lim [ d 2 r~?(r)Vm>(r) (20) 

R— *0 J 

where the regularized forms defined below are to be employed for 5(r) and \&(r) as 
well as for Q. 

We solve the eigenvalue problem of the Hamiltonian 

(E',a',m'\H a - E\E,a,m) = 0. (21) 
The solutions have to be normalizable, so we impose 

hm| d 2 r^ m {r)^ E ^ m {r) = 1. (22) 

For a solution of ( pip we use the continuous ansatz ( |T6| ) as a regularized form. 
Putting this in ([H]) we get a contribution only from the discontinuity of the first 
derivative at f = R, which has to vanish for R — > 0, such that 



lim 

R^O 



f=R 



f=RJ 



0, (23) 



ensuring hermiticity of H a for R — > 0. This condition is, of course, fulfilled by 
demanding differentiability for all finite R, which is what the second approach in 
section 3 was based on. Allowing, however, the expression in the square brackets 
in ( p3|) to be finite for any R > and only demanding that it vanish in the limit 
R — > is a weaker condition on the eigenfunctions. 

A detailed investigation of ( |2"3"| ) using the ansatz QT5D may be summarized as 
follows. There are essentially two ways for an eigenfunction to fulfill (|23|). The 
first is that the leading power in R is great enough to make both terms within the 
square brackets in (^) vanish separately as R —>■ 0. The second is that leading 
non-vanishing contributions cancel each other, while higher order terms still vanish 
separately. This leads again to the condition ( |IBD to be fulfilled. Since it cannot 
be fulfilled for both partners of a supersymmetric pair simultaneously, it can be 
fulfilled either for a supersinglet, i. e. an eigenfunction for E = 0, or for a partner 
of an eigenfunction that fulfills (|23|) in the first way. 

Besides normalizability and ( |2"3"D we require E > and supersymmetric pairing 
of the eigenfunctions for E > 0. 

The pairing implies that an eigenfunction for E > fulfilling fl23"| ) is admissible 
only if its partner also fulfills (|23|). But we have to take into account that the su- 
peroperators contain a derivative in f and will in general map a continuous function 
with a discontinuous first derivative at R to a function that is discontinuous at R. It 
is, however, always possible to restore continuity by adding a contribution that does 
not change the norm, i. e. a representative of the null vector. Moreover, the behavior 
for r < R can always be arranged to have the eigenfunction in the form (0), which 
we will call the basic form of the eigenfunction. The freedom of representing the null 
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vector is an important feature of the regularization and will be discussed in more 
detail in the next section. 

The normalized eigenf unctions of Laguerre type are 



Aa/7T 



(E-a- 1/2 -m- a)! 
T(E- a +1/2) 



6(f - R)r m+a e^Ltt%2^ a (r 2 ) 



+6{R - r) 



L Fi(m9{m) + a + 1/2 



R 2 E 



1 + \m\,R 2 + 



a: 



xe -(i+ a /R 2 )f 2 /2 f H iFi ( m £( m ) + a + 1/2 



R 2 E 

R 2 + a' 



a 



1 + \m , 1 + f 
i? 2 



e imip (a 
(24) 

for i£ — er — 1/2 — m — a = 0, 1, 2, ... with m + a > for o = — 1/2 and m < 0, and 
m + a > —1 for o = +1/2 and m > 0. 



(£-<7-l/2)! 



Av^F\ T(E - (7 + 1/2 - m - a) 



6(r - R)f 



~ — ni — a 



e ^E-a-i/zv ) 



xe 



+9(R - f) 



-{l+a/R 2 )r 2 /2~\m\ 



R- m - a -\ m \e^L { E m a _^ /2 {R 



i F x (md{m) + a+ 1/2 



f |m| lFl(m0(m) + <7 + 1/2 



i? 2 +C 
£ 2 £ 

i? 2 + a' 



1 + |m|,i? 2 + a) 

L + Irak [1 + S- ) r 2 
1 ''V i? 2 / 



C<7 



(25) 

for _E - a - 1/2 = 0,1,2,... with m + a < for a = +1/2 and m > 0, and 
m + a < 1 for a = — 1/2 and m < 0. The eigenf unctions in fl24]) and (|25|) contain 
all eigenfunctions of the Hamiltonian by the direct approach. But (^) still allows 
further eigenfunctions of Laguerre type along with their superpartners. This reduces 
the vacancy in the spectra by one (m+cr)-value in comparison to the direct approach. 

The other eigenfunctions fulfilling (|23|) are supersinglets. In order to be nor- 
malized as R — > they acquire an independent normalization factor that lowers 
the leading power in R in (f23|). This very fact ensures that fulfilling ( p3|) in lowest 
order is sufficient for (^) to be fulfilled to all orders, since the higher powers vanish 
separately. We call the eigenfunctions with an independent normalization factor 
singular eigenfunctions. 

For E = and a = —1/2 we find the normalized singular eigenfunctions for 
m + a = 1 



E=0,a=-l/2,m 



1 



A 



-7T 



lni? 2 



x 



9(r - R)f- X e^- + 9(R - f )R- 1 e a/2 e- {1+ ^ ) ^ ( ^ 

\R 



C-l/2 



(26) 
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and for m + a > 1 

*ko i(r =-l/2,m( r » 



m+Q— 1 



A^tt(^ t + a— i e -[r(-m + 1) - T(-m + 1, a)}) 



x 



- R)r m - a e^- + 6{R - f)R- m - a e a/2 e- {1+ ^ ] 



a_i£i ( r 

2 



1? 



e^C-1/2 (27) 



with the incomplete gamma function We have thus recovered all the correct 
eigenfunctions found within the approach in 3.2. 

For E = and a = +1/2 we find still another type of singular eigenfunction. 
For m + a = — 1 we obtain 



E=0,(T=+l/2 



,m(r,</>) 



1 



A v — 7T In i? 2 



x 



- J R)f" 1 e^r(l, f 2 ) + - f) J R- 1 e^r(l, R 2 ) e - a/2 e {1+ ^ (4 

\R 



and for m + a < — 1 

*E=0,<x=+l/2,m( f > ¥>) 



e im ^C+i/2 
(28) 



-m— a— 1 



A xMz^zt + e- a (-«)-" l - 1 [r(m + 1) - T(m + 1, -a)]) 



(f - R)f m+a e 



f£ r(— m — a, r 
r(— m — a) 



r.2\ 



+6(R - r)R m+a e^ 



r 2 r(-m - a,R 2 



Tf— m — a) 



. e -«/2 e (l+^)V 



1? 



e imv C + i/2. (29) 



They were not anticipated in section 3 by either approach. 
The spectrum of H a is displayed in figure 1. 



4.2 The Hilbert space in the regularized coordinate repre- 
sentation 

The scalar product in the regularized coordinate representation is given by 

(E|tf) = lim f d 2 rZ*(r)^(r), (30) 

R — >0 J 

exhibiting the necessary properties. An important difference between the usual 
coordinate representation and the regularized coordinate representation is that in 
the former the null vector is represented by the function that vanishes everywhere, 
while in the latter also functions N(r) that fulfill 

lim f d 2 rN*(r)N(r) = (31) 

R— >0 J 
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also represent the null vector. We restrict the set of functions representing the null 
vector to those that fulfill lim^o N(r) = excluding probability densities that 
are different from zero on curves and points in the plane. Along with the class of 
functions that represent the null vector, every state is represented by an equivalence 
class of functions by addition of representatives of the null vector. The corresponding 
probability density defined by lim^ ^*{r)^{r) is unique by using the basic form 
of a wave function defined below. 

The linear space of functions with the above scalar product is still not a Hilbert 
space, since we have to ensure completeness; i. e. every Cauchy sequence has to 
converge. To this end we specify the wave functions \P(r) G 7i a - The numerable set 
of eigenfunctions of H a (]2~4])-(|2~9|) constitutes an orthogonal normalized system with 
respect to the scalar product 

(E',a\m'\E,a,m) = \^mj d\^ E , alm ,{r)^ E ^ m {r) = 5 E i E 5 a ' a 5 m > m . (32) 

We obtain general wave functions from linear combinations of the eigenfunctions 

*W = C E,a,m^E,a,m(r) with C E ^ m C E ,a,m < OO. (33) 

E,a,m E,cr,m 

This defines the basic form of a wave function as the sum of eigenfunctions in the 
basic form. Wave functions will always be given in the basic form, unless otherwise 
indicated. The space of functions obtained in this way is a Hilbert space by con- 
struction, since we define its elements as limiting objects of the Cauchy sequence of 
an increasing number of members of the sum in ([33]). 

The basic form of a wave function Sb as ic(r) * s obtained from an arbitrary equiv- 
alent wave function 2 non _ basic (r) by 




This relation also expresses the fact that the eigenfunctions of the Hamiltonian 
supply a resolution of unity. It is this requirement that enforces the construction of 
the Hilbert space as presented here. 

For a = the Hilbert space is clearly L 2 <8> (C 2 . Also for non-integer a ^ the 
eigenfunctions of Laguerre type provide a resolution of unity for L 2 <S> C 2 P3"fl . For 
| a | > 1 the singular eigenfunctions are added. For integer a^O from the functions 
of Laguerre type those with (m + a) of the vacancy are removed. 

The projection of wave functions CE,a,m = (E, a, m\ty) in the regularized coordi- 
nate representation on the elements of the complete orthogonal system yields a ma- 
trix representation where each state is uniquely represented by the vector ({cE )(J)m }), 
e. g. the null vector is given by (0, 0, ...). 

The pairing of states by supersymmetry is possible in the regularized coordinate 
representation only by the freedom of adding a representative of the null vector. 
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E. g. the supersymmetric pairing of ^ E,a=-i/2,m(f) with E,a=+i/2,m-i(f) ma Y be 
expressed by 

lim / d 2 r~*(r)Q a y Eta= ^ 1/2tm (f) = lim / d 2 r H*(f)[Q A iff= _ 1/2jm (f) + N(f)} 

= lim / d 2 rS*(f)(±V^)* £)(T = + i/2 lTO - 1 (f), (35) 

with the upper sign for eigenf unctions from ( ^4f) and the lower sign for eigenfunctions 
from (Po]). The function Q Q ^/_B,o-=-i/2,m(^*) with a step at r = i? is repaired in order 
to obtain ^ E,a=+i/2,m-i{'r) m the basic form by addition of the appropriate N(r). 



4.3 Probability densities 

At this point it is useful to discuss the modulus squared \E'*(r)\E'(r*) in the limit 
R — > related to measurement as a probability density for the eigenfunctions of the 
Hamiltonian. 

Of course, for the eigenfunctions of Laguerre type the functional form of the 
exterior part of the modulus squared extends to all r ^ if R — > 0. So only the 
single point r = remains to be considered. 

For | a | > 1 we find that the probability density is zero at r = 0. The vast 
majority of the eigenfunctions of Laguerre type vanish as r — > 0, yielding a density 
that is continuous everywhere. We have a non-continuous behavior at r = if the 
outside density reaches a constant value or diverges as r — > 0, since the probability 
density is still zero at r = 0. Thus, for |a| > 1 the probability densities of the 
non-singular wavefunctions built on eigenfunctions of the Laguerre type avoid the 
location of the flux tube. 

By contrast, for \a\ < 1 the value at r = agrees with the behavior of the exterior 
part as r — > 0, such that either the probability density is continuous everywhere or 
the value at r = diverges as R — > along with the exterior part diverging as 
r -»• 0. 

The normalization of the singular eigenfunctions implies 

lim / d 2 r^ m (r)^ m (r) = 1 (36) 

for any Q containing r = 0. If Q did not contain r = 0, the integrand would converge 
uniformly to zero and the integral would vanish as R — > 0. This means that the 
modulus squared of a singular eigenfunction for finite R is a delta convergent function 
|23, i.e. 

hm n=0,a,m(^E=0,a,m(^ = (37) 

where the limit R — > is to be taken outside of an integral. Thus, the probabil- 
ity density of a singular eigenfunction is that of a two dimensional point particle, 
concentrated at the location of the flux tube. 

Singular eigenfunctions occur for |a| > 1 just when the probability densities 
of non-singular wavefunctions avoid the flux tube as mentioned above. Although 
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not being formally the reason, the fact that the probability density of the singular 
eigenfunctions and the eigenfunctions of Laguerre type are attributed to disjoint 
regions in the plane conforms nicely with the mutual orthogonality. Indeed the sin- 
gular eigenfunctions are orthogonal to any function H(f) yielding a Taylor expansion 
about r = 

hm J d 2 r *£ =0)CT , m (r)E(r) = 0. (38) 

From this we also conclude that ^ s E=Qam {r) itself neither is a delta convergent func- 
tion nor does it converge to a derivative of the delta function. Since this exploits all 
possibilities for generalized functions concentrated at a point ||21|| , the corresponding 
functional must be zero. 

We define a general singular wave function \E ,s (r) as a linear combination of 
singular eigenfunctions. It also has the properties (|37|) and ([38]) if it is normalized. 



4.4 Classification of the Hilbert spaces 7i a 

By a a sequence of Hilbert spaces is parametrized. We have determined each Hilbert 
space by solving the eigenvalue problem of the Hamiltonian exactly for all values of 
a. It is, therefore, not necessary to resort to perturbation theory. Still, in view of 
extensions of the system and analogous situations the question arises whether we 
can get arbitrarily close to the exact solution of a system with a' = a + Aa by 
perturbation theory in Aa if the system with a is exactly known. 

A necessary condition for this is that exact solutions of the perturbed system 
can be expanded in terms of the exact solutions of the unperturbed system. An 
exact solution of the perturbed system that is orthogonal to all exact solutions of 
the unperturbed system can never be reached by a perturbation series. 

The real numbers that a can assume decompose into the isolated points ±1; ±2; ... 
and the complementary open intervals (—1,1); (1,2), (2,3),...; (—2,-1), (—3,-2),... 
This is a result of the occurrence of the singular states and the vacancy in the 
(m + a)-range for a = ±1, ±2, ... These sets represent equivalence classes with re- 
spect to the relation TC a = TC a '. For any a and a' in different equivalence classes we 
have H a ^ H a ' such that there is at least one state in TC a or TC a ' that is orthogonal 
to all states of the other. For any two neighboring equivalence classes it turns out 
that both Hilbert spaces contain at least one state that is orthogonal to all states 
of the other. For these states perturbation theory fails and the above equivalence 
classes constitute perturbatively disjoint sectors. We summarize: 

i) As expected, Aa may not be arbitrarily large for a successful application of 
perturbation theory; in fact |Aa| < 2. 

ii) However, very small Aa does not guarantee successful application of perturbation 
theory, since a and a + Aa could still be in different equivalence classes. 

iii) The open interval (—1, 1) is distinguished as the interval length is exceptionally 
2, while for the other equivalence classes we have interval length of 1 or 0. a = 
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is the only integer that is useful as a perturbative starting point for its neighboring 
values. The Hilbert space is that of square-integrable 2-spinors 



The dependence of H a and the superoperators on a does not anticipate the de- 
composition into disjoint perturbative sectors. However, by a we parametrize mag- 
netic flux, which is the global topological invariant associated with a U(l) principal 
fiber bundle characterizing the topology of the bundle Thus, the discontinuous 
behavior of the Hilbert spaces reflects the differing topologies, i. e. discontinuous 
geometries, as a is varied. 

4.5 An index for the singular flux tube 

Index theorems relate the spectrum of a differential operator on a manifold to a 
global topological invariant. The number of eigenstates of E = with a = —1/2 



minus those with a = +1/2 yield an index, referred to as the Witten index |f20| [25 
in super symmetric systems. It corresponds to the index of the euclidian Dirac 
operator in two dimensions, which is subject to index theorems. The Atiyah-Singer 
index theorem || for manifolds without boundary and the Atiyah-Patodi-Singer 
|TC)|j index theorem for manifolds with boundary apply to compact manifolds. For 
the U(l) case in two dimensions the Atiyah-Singer index theorem reduces to the 
statement that the above index equals the total flux in units of the elementary flux 
quantum. The index theorems are not applicable to our system as the plane is not 
a compact manifold. Both the above index and the magnetic flux through the plane 
are infinite due to the presence of the homogeneous magnetic field. 

More can be said in view of the fact that in the limit R — > both the eigenfunc- 
tions of E = and the flux may be attributed to disjoint subspaces of the plane: 
the point r = and its complement. 

For r ^ the eigenfunctions of Laguerre type with E = and o = —1/2 that 
are not localized close to the flux tube can be attributed to flux quanta of the 
homogeneous field as for a = [0. 

For r = we define the index I s (a) as the number of singular eigenfunctions 
with a = —1/2 minus the number of singular eigenfunctions with a = +1/2. If [a] 
denotes the next integer below a, it is given by 

I s (a) = [a]0(a) - [-a]0(-a), (39) 

displayed in figure 2. For integer a we have I s (a) = a, as if we had a compact 
manifold as in [51. 



5 Explicit breaking of supersymmetry by g ^ 2 

Concerning the applicability of the results obtained so far, we are facing the fact 
that already in vacuo g ~ 2.0023 by QED corrections. Even by this tiny deviation 
from g = 2 supersymmetry is explicitly broken. The question is to what extent our 
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results remain valid for g ^ 2. We discuss only the effect of QED corrections here, 
but for other perturbative background effects the arguments might be analogous. 

The deviation of g from 2 is calculated under the assumption of a homogeneous 
field that is sufficiently weak, i. e. | - jl ■ B\ < Mc 2 p7fl[pg|. For the outside 
region we can assume this to be true. By contrast, inside of the flux tube the 
magnetic field grows as R — > 0. Assuming that R — > idealizes R ~ ja^Ac, we 
have \ —jl-B\ « Mc 2 and it is not justified to put g ^ 2 there. The most reasonable 
choice is to maintain g = 2 within the flux tube. This may also be rephrased 
heuristically following the correspondence principle. A strong magnetic field implies 
high occupation numbers of the quantum electromagnetic field. The limit of high 
occupation numbers corresponds to classical behavior. As a result, contributions 
to QED processes from the fluctuations of the quantum electromagnetic field are 
suppressed. 

The most obvious implementation of g = 2 inside and g ^ 2 outside would be to 
couple the entire magnetic field with g = 2 for r < R and with g ^ 2 for r > R. The 
result would, however, be inconsistent with our regularization. On the one hand, for 
| a | > 1 the expectation values of the Hamiltonian with the singular eigenfunctions 
would not reproduce the eigenvalue. On the other hand, the regularized Hamiltonian 
would not yield the Hamiltonian of the homogeneous field for a = 0. Therefore, we 
proceed differently. 

We couple Be z by g = 2(1 + k) ^ 2 everywhere and (a§/irR 2 )e z by g = 2 
inside of the flux tube, which amounts to a term (1 + k)BS z + 9(R — r)(a^/iiR 2 )S z 
in the Hamiltonian. aQ/nR 2 dominates B in the limit R — > and we effectively 
attribute to the interior g = 2 and to the exterior g ^ 2 in a consistent way. The 
new Hamiltonian is 

H« = H a + kS z = {Q a , Qi} + kS z , (40) 

entailing the explicit breaking of supersymmetry 

m,Q a ] = KQ a and [H^QI] = -kQI (41) 

A lower bound for the expectation values of H£ is now {^/\H^\^/) > —k/2 in- 
stead of zero. For the solution of the eigenvalue problem of H* we have still 
[H*, L z ] = [H£, S z ] = 0. The matching condition (|23| ) is not affected, since the 
leading contributions at small R of both the outside and the inside parts of ([16]) are 
independent of k. We obtain the eigenfunctions for k ^ from those for k = given 
by (]24])-(]29D replacing E by E K — kg. Thus as functions of r*they are unchanged, 
while the eigenvalues E K = E + na of H* are shifted. 

The eigenfunctions with E K = kg are now the supersinglets. The other eigen- 
functions are still paired by Q a and Q^, although having different energy eigenval- 
ues. One might object that, since the superoperators no longer commute with the 
Hamiltonian, the pairing is not enforced as for the k = case, where the pairing has 
effectively determined the correct "boundary conditions" . Then additional unpaired 
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eigenfunctions could enter the spectrum for k^O. This is not true, however. Ad- 
ditional eigenfunctions would have to disappear abruptly in the limit k — > 0. Since 
k 7^ is a perturbative effect, one has to impose a continuous limit k — > 0. 



6 The singular Aharonov-Bohm system 

The special value of B = has not been covered so far, as for r — > oo we can 
no longer impose the vanishing of the eigenfunctions of the Hamiltonian. By the 
discussion of the last section only the case of g = 2 needs to be considered here. 

Lacking a knowledge of the behavior of the eigenfunctions as r — ► oo, we have 
to start with a more fundamental criterion for selecting the correct eigenfunctions. 
Ultimately, the set of correct eigenfunctions of the Hamiltonian has to supply a 
resolution of unity in the Hilbert space. 

We consider the special case a = 0, where this system reduces to the free particle 
in two dimensions. The corresponding eigenfunctions are plane waves, not actually 
constituting a normalized orthogonal set. Still these eigenfunctions supply a res- 
olution of unity in the Hilbert space of square-integrable 2-spinors 7io = L2 £S> (D 2 
by theorems on the Fourier transform. Obviously, the behavior of the system for 
r — > 00 does not affect the Hilbert space. A normalizable orthogonal system span- 
ning 7io = L2 ® (D 2 is given by the eigenfunctions of the Hamiltonian of the homo- 
geneous magnetic field. 

This can be generalized to a ^ 0. We demand that the solutions of the eigenvalue 
problem of H^ B at B — yield a resolution of unity in the Hilbert space Tt a spanned 
by the normalized orthogonal set of eigenfunctions (p4|)-(p9|) at B ^ 0. For a / 
the system looks the same as for a = as r — > 00. Generally, the value of B by 
which the systems might differ as r — > 00 does not affect the Hilbert space. But we 
have shown that the Hilbert space underlying the system does depend on the value 
of a, i. e. on the behavior at the origin. 

According to 4.4 any system with B 7^ and a' within the same equivalence class 
as a would supply a suitable reference basis. We choose a = a' for convenience. 

The operators 



4 \r dr dr r 2 dip 

+ 0(r~R) f-Sll + ^h ^-r)^[^r + — , S- ) . (12) 





2r 2 dip 

(O* = S-(0* = S-^- ~ - ~ + - + a9(R - r) (± - i 

dr r dip r \K Z r 



(43) 



are obtained by multiplication of (0) and (|i"5|) with A 2 and A l , respectively, and 
the subsequent limit A — > 00. These operators are now measured in units of 2h 2 /M 
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and 2 1 / 2 h/(M) 1 / 2 , respectively. The supersymmetry algebra (Q) is fulfilled by Q^ B , 

(Q : y ! y and n,y>. 

We solve the eigenvalue problem of the Hamiltonian in the form 



k 2 



(E,a,m\ #r- T j\k) = 



(44) 



A necessary condition on the solutions to yield a resolution of unity is that they have 
to be within the Hilbert space; i. e., all (E, a, m\k) must be finite and at least one for 
a given \k) must be different from zero. Hermiticity of the Hamiltonian imposes a 
further condition on the solutions, which requires orthogonality (k'\k) = if k 7^ k! . 
Due to [H* B } L z ) = [H£ B , S z ] = 0, we use the ansatz * fci<T>ro (r, <p) = W,r»e^C- 
Inside of the flux tube we have again to solve a Kummer equation. Outside of 
the flux tube we have a Bessel equation for k > 0, which reduces to the Laplace 
equation for k = 0. If we impose continuity at r — R, from the eigenvalue equation 
(fyj) the matching condition ( f23|) follows again. We will not combine the two linear 
independent solutions for k > to yield different iable solutions at r — R, such 
that (p3|) would be fulfilled even before taking the limit R — > 0. Instead, in order 
to proceed coherently with the case of the additional homogeneous field we only 
demand continuity. The discussion based on the leading behavior of the solutions in 
section 4.1 is repeated literally by using only one of the linear independent solutions 
for r > R. The result is for k > 



k,a,m 



{r-R)J m+a {kr)+6{R-r) 



R-\ m \e%J m+a (kR) 



1 F 1 (m9(m) +a + l/2- « 1 + \m\, a 



x 



e- ar2/2jR2 r |m| iiq(m0(m) + a + 1/2 - 



2JU2 



R 2 k 
4a 



1 , a 



tmip 



(45) 



with m + a > for a 



k,a,m 



-1/2 and m < 0, and m + a > —1 for a = +1/2 and m > 0; 

R-^e^J^ a (kR) 



{r-R)J_ m _ a (kr)+6(R-r) 



l F L (m0(m) + <7+l/2- 



4a ' 



21.2 



X 



e-^^r^ljF^m^m) + a + 1/2 - 1 + |m|, -^-r 2 ) 



4a 



i? 2 



D mup 



1 + Iml, a) 



(46) 



with m + a < for a = +1/2 and m > 0, and m + a < 1 for cr = — 1/2 and m < 0. 
For = the solutions are 



fc=0,CT==Fl/2,m 



i/[0(r - R)r 



q=(m+a) 



+ 0(i2 - r)R Ta e ±a/2 e Tar2/2R2 r Tm ]e imv C Tl/2 



(47) 



for a = +1/2, ±m < and +(m + a) > 1/2. There are three types of solutions for 
k = concerning their normalization behavior. 
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i) We have v = ij±(™+«H^V2|( ± ( m + a )_ 1 )-i + ( ±a )±™-i e ±« 7 ( Tm + i )±a )]-i/2 
for + (m + a) > 1 and o = +1/2. The normalized eigenf unctions can be found by 
directly performing the limit of the vanishing homogeneous field, i. e. A — > oo in the 
normalized eigenfunctions in and (|29|) . The modulus squared is still a delta 
convergent function. 

ii) We have v = [— n lni? 2 ] -1 / 2 for +(m + a) = 1 and a = + 1/2. The eigenfunctions 
occur for integer a^O and they correspond to the only eigenfunction in the (m + a)- 
vacancy in the case of the additional homogeneous magnetic field. However, they 
are not obtained by the limit A — > oo from the normalized functions (|26|) and (p8[) , 
since they yield a logarithmically divergent norm, v is determined by demanding 
that {^e=o am\^k=o~am) = 1 f° r arbitrary A ^ oo. Still as for finite A, the modulus 
squared yields a delta convergent function, since the weak logarithmic divergent 
behavior as r — > oo is compensated by any decaying test function. 

iii) For 1 > + (m + a) > 1/2 and o = +1/2 we put v = 1 for convenience. These 
solutions share their property of not being normalizable with the continuum states 
with k > 0. 

The eigenfunctions P5D-(pE7D constitute a resolution of unity in Ti a . This follows, 
on the one hand, from the fact that the eigenfunctions of the free particle represent 
a resolution of unity in L 2 <S> (D 2 . In the case of integer a ^ both for finite A and 
A —>■ oo the vacancy occurs for the same (m + <j)-value. On the other hand, the 
singular eigenfunctions for A — > oo occur for the same quantum numbers m and a 
as for finite A. They yield (^ r */ A 7n 7 00 |^ r J^) = 8m,m'5ay for all finite A and, therefore, 
in the finite subspace of the singular wave functions we have a resolution of unity. 
Since the singular eigenfunctions and the non-singular ones are mutually orthogonal 
by fl38j) , the corresponding resolutions of unity do not interfere and the sum of both 
yields a resolution of unity of the entire TC a . 

Having only a finite number of eigenfunctions with k = 0, we define I (a) as 
the number of eigenfunctions for k = and a = —1/2 minus those for k = and 
a = +1/2. It is given by 

in for a <E [n - 1/2,71 + 1/2) n = -l,-2,... 
I AB (a) = J for a G [-1/2, 1/2] (48) 
{ n for a e (n- 1/2, n+ 1/2] n = l,2, ... 

sketched in figure 2 on the right hand side. -P(a) is the same as for the nonvanish- 
ing homogeneous field. (|48|) coincides with the Atiyah-Patodi-Singer index theorem 
for compact manifolds with boundary, according to Jl4| by use of the appropriate 
boundary conditions. For integer a both I s (a) and I AB (a) give the same value 
a in agreement with the Atiyah-Singer index theorem || valid on compact mani- 
folds without boundaries, which is plausible since in this case a compactification is 
possible. 
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7 Concluding remarks 



Our results do not depend on our special arrangement of the magnetic field, i. e. 
regularization. Smearing out the discontinuity of the magnetic field between R — 77 
and R adds to the matching condition (|23| ) a correction of 0{rj) that is forced to 
vanish by the limit R — > 0. 

The need for a regularization is dictated by consistency. And for this reason also 
the consequences of the use of a regularization are to be taken seriously. Ultimately, 
in what way these idealized structures could be realized in nature is yet another 
question to be resolved by experiment. 

In view of the quantum Hall effect for a = 1, 2, ... the eigenfunctions of Laguerre 
type preserve the structure of the lowest Landau level of the homogeneous field 
a = supplying the states to be occupied for the incompressible quantum fluid. 
The singular states are orthogonal to these, thus not participating. If the singular 
states remained unoccupied, the number of electrons per non-singular flux quantum 
would be the same as for a = 0. However, the total magnetic flux would be increased 
by the contribution of the singular flux tube. In other words if there were a singular 
flux tube without occupation of the corresponding singular states, the filling fraction 
would be effectively decreased. 

Finally, this system might be interesting as a toy model for supersymmetric field 
theories, for two reasons. First, because the singular flux parametrizes different 
topological situations and thereby controls the subspace of supersinglets, related 
to spontaneous breaking of supersymmetry. Second, the deviation from g = 2 
provides an explicit breaking of supersymmetry from effective physics such as QED 
corrections at lower scales Ac than the scale A of the quantum mechanical system. 
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FIG. 1 (previous page). Spectrum of H a 

A state of Laguerre type with a = — 1/2 is indicated by -,- and a state of Laguerre type 
with cr = +l/2by j-. A supersymmetric pair of such states is indicated by +. Singular 
states with a = —1/2 are indicated by u and singular states with a = +1/2 by n. For 
a = the eigenvalues of H a constitute a regular rectangular lattice within a 3/8-sector 
of the E-(m + <r)-plane. The rectangular lattice acquires a line defect for a / 0. For 
integer a / the defect is a vacancy line between the left and the right part of the 
spectrum, i. e. one (m + cr)-value is missing, entailing the absence of the corresponding 
Fourier components ~ e imi f among the eigenfunctions of the Hamiltonian. For non-integer 
a there is a vertical shift of the right block against the left block by the non-integer part 
of a. All non-integer values of a between two integers show the same qualitative pattern 
differing only in the amount of vertical shift of the right block. The spectrum is given for 
half- integer values. 
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FIG. 2. The dot • indicates the value of 1(a) at the steps. 
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